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The need for strategies able to accurately manipulate quantum dynamics is ubiquitous in quan- 
tum control and quantum information processing. We investigate two scenarios where randomized 
dynamical decoupling techniques become more advantageous with respect to standard deterministic 
methods in switching off unwanted dynamical evolution in a closed quantum system: when deal- 
ing with decoupling cycles which involve a large number of control actions and/or when seeking 
long-time quantum information storage. Highly effective hybrid decoupling schemes, which combine 
deterministic and stochastic features are discussed, as well as the benefits of sequentially imple- 
menting a concatenated method, applied at short times, followed by a hybrid protocol, employed 
at longer times. A quantum register consisting of a chain of spin-1/2 particles interacting via the 
Heisenberg interaction is used as a model for the analysis throughout. 



I. INTRODUCTION 

The constructive role of randomness in physical processes has been demonstrated in various areas of research. In 
stochastic resonance |lj, f° r instance, a weak signal can be amplified by the assistance of an appropriate noise. In 
quantum information processing, noise can intensify the speed-up of quantum walks over classical ones 0, dissipa- 
tion may offer new possibilities to implement gate operations in quantum computing Q, while static perturbations 
characterizing faulty gates can enhance the stability of quantum algorithms In quantum communication, the use 
of random operations decreases the communication cost of achieving remote state preparation and of constructing 
efficient quantum data- hiding schemes j>j . Finally, random unitary operators have been recently suggested as allowing 
efficient parameter estimation for open quantum systems . 

In the context of coherent quantum control, the advantages of stochasticity have only recently been addressed 
0, H El El Oi • Within the framework of dynamical decoupling methods, in particular, analytical bounds derived 
in pointed to situations where randomized protocols are expected to outperform their deterministic counterparts 
in suppressing unwanted unitary dynamics as well as decoherence in open quantum systems. The idea of merging 
together deterministic and randomized designs into hybrid control schemes, where benefits from both approaches may 
be simultaneously exploited, was also proposed in in general control-theoretic terms, and independently validated 
in illustrative situations in (see also (l^l- 

Here, we focus on exploring the advantages of randomization in establishing efficient control schemes for arbitrary 
quantum state stabilization, that is, for engineering a quantum memory Efficiency is assessed in terms of both the 
number of control operations needed to achieve a desired fidelity level and the rate at which residual errors build 
up in the long run. We show that by interpolating the most effective (deterministic) scheme known for short-time 
decoupling with the best available (randomized) scheme for long time, very high performance over the entire time 
axis may be ensured. 

II. SYSTEM AND CONTROL SETTING 

Dynamical decoupling (DD) techniques have been extensively discussed both in the original high-resolution nuclear 
magnetic resonance (NMR) setting and, more recently, in connection with robust quantum information processing, 
see for instance 0, 0, El f° r representative contributions. Assume, for simplicity, an isolated finite dimensional 
target system S described by a (possibly) time dependent Hamiltonian Ho(t). The basic idea of DD is to modify the 
dynamics of S by adding to Ho(t) an appropriate time-dependent control field H c (t). The overall propagator under the 
total Hamiltonian and the control propagator in units h — 1 are, respectively, U{t) — Texp[— i f Q (Ho(u) + H c (u))du] 

and U c (t) = Tcxp[— i J Q H c (u)du], where T indicates time ordering. A transformation to a logical frame that removes 
H c {t) is commonly performed, leading to a controlled evolution described by 
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il{t) = Ut{t)U{t)=T 



exp 



i I Ho(u)du 



(1) 



where Ho(t) = U£(t)Ho(t)U c (t) is the logical Hamiltonian. When Ho(t) is time-independent and the perturbation 
is cyclic with cycle time T c , i.e. U C (T C ) = i, where 1 is the identity operator, physical and logical frame coincide 
stroboscopically at T n = nT c , n £ N. Using the formalism of average Hamiltonian theory (AHT), the evolution 
operator in the logical frame may be expressed as U(nT c ) = exp(— iHonT c ), where 



U{T C ) =Texp 



Ho(u)du 



- e- iHoT < (2) 



defines the average Hamiltonian Hq — J^kLo -^o*^> with each term computed from the Magnus expansion (l7lll8|. 
A sufficient convergence criterion for the series is given by kT c < 1, where k — \\Hq\\2 and ||^4||2 = max |eig(.A)|, 
WA = A*. 

The above time average for Hq may be conveniently mapped into a group-theoretic average. In the framework of 
bang-bang DD, in particular, control actions correspond to arbitrarily strong and effectively instantaneous rotations 
successively drawn from a (projective representation of a) group, Q = {gj}, j = 0, . . . , \Q\ — 1. The propagator at T c 
is written as 

ISI-i 

U(T C ) = U(T C ) = Y[ g\U {t k+1 ,t k )9k , (3) 

fc=0 

which translates into a cyclic sequence of pulses = fl'fcpt-p k = 1, . . . , \Q\, separated by intervals At = tk+i — t% of 
free evolutions, leading in turn to a cycle time T c — \Q\At. The zeroth order contribution of the Magnus expansion, 

which dominates in the limit T c — > 0, is therefore = Hg = \G\~ X J2j sJ-^oSj- 

The so-called time suspension is the DD goal we focus on here, that is, we want to develop pulse sequences able 
to approximate the evolution operator as close as possible to 1. How well we succeed at preserving a given initial 
state is reflected, for instance, in the proximity of the input-output fidelity F(T) to its maximum value 1, where 
F m (T) = MU(TM\\ 

A deterministic protocol based on a fixed control path of a representation of Q and aiming at achieving first-order 
decoupling, — 0, will be referred to as "periodic deterministic decoupling" (PDD) protocol We will assume 
here that the first PDD pulse occurs only at t\ — At, that is, go = i. Its simplest stochastic version is obtained 
by randomly picking elements over Q, such that the control action at each t n = nAt [to included] corresponds to 

p(r) = 

9i9p i,3 = 5 . .. , \G\ — 1- This leads to what we call "naive random decoupling" (NRD) - an intrinsically acyclic 
method. Bounds on the worst-case pure-state expected fidelity at time T, F(T) = min^\ F^(T), were established 
in Ref. 0. For PDD, in the limit TT c k 2 < 1, we have: F(T) > 1 - C(T 2 T c 2 k 4 ), while for NRD and TAtn 2 < 1: 
E{F(T)} > 1 — 0(TAtn 2 ), where E denotes ensemble expectation over all control realizations. We note that the 
bound for NRD still holds in the case of a time dependent Hamiltonian as far as 1 1 i^o (^) 1 1 2 is uniformly bounded 
in time by k > 0. Within their regime of validity, these bounds indicate that NRD should outperform PDD when 
\Q\ 2 (T AtK 2 ) 3> 1, which is the case when large control groups and/or long interaction times are involved. 

In practice, we avoid the extremization procedure required to determine the worst-case pure-state expected fidelity. 
Instead, in order to get a state-independent estimate of the performance of different DD methods, we invoke gate 
entanglement fidelity, F F (T) poj ]. This quantity is linearly related to the average input-output fidelity over all pure 
states of the system [2l|, |23], and it may be computed as F e (T) — |Trace(J7(T))/d| 2 , where d is the dimension of 
the system's state space. Our control objective is then to get as close as possible to E{F e (T)} — * 1. In the Monte 
Carlo simulations to be presented below, ensemble expectation, E{F e (T)}, is further replaced by an average over a 
sufficiently large statistical sample of control realizations, leading to what we denote ((F e )). 

In order to concretely illustrate the benefits of randomization, we concentrate on a relatively simple, yet physically 
relevant, example - that is, to completely refocus the internal evolution of a chain consisting of N strongly coupled 
spin- 1/2 particles (qubits) described by the Heisenberg model, 

N N-l 

Ho = Y,\ 1 + J Y, ( X i X i+i + Y ^+i + A^^+i) ■ (4) 
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Here X, Y, and Z denote Pauli operators, u>i is the frequency of qubit i, J is the coupling parameter, and A determines 
the anisotropy. Only nearest-neighbor interactions are considered, which is a fairly good approximation for couplings 
exponentially decaying with the qubit distance - as arising, for instance, in quantum dot arra ys 23l| - or decaying 
cubicall y - as it is the case for dipolar interactions of NMR crystals and liquid-crystals 0, l24j. or electrons on 
Helium |25|. 

We consider qubits with approximately the same frequency a->; ~ lu. Accordingly, in order to remove the phase 
evolution due to the one-body Zeeman terms, we perform a transformation to a frame rotating with frequency 
uj and characterized by the operator Unit) — exp 

H w *Ef Zi/2]. We thus work in a combined logical-rotating 
frame, whereby the effective Hamiltonian becomes Hn(t) = [/j(f)/7j ? (t)[77o(t) — Zi/2]Un(t)U c (t). While this 
approximation is accurate for a class of physical systems (notably, homonuclear NMR samples), the restriction is not 
fundamental. If the spread of the single-qubit is significant (so that a common rotating frame does not exist), schemes 
capable of additionally refocusing the Zeeman terms may be constructed without adding to the overall complexity of 



Following the general idea of |15l j26j, a PDD protocol capable of refocusing the nearest-neighbor couplings of 
Eq. @ for arbitrary parameter values may be built by recursively nesting DD sequences based on the group Qi = 
{li, Zj, Xj, Yi] for each even qubit, i = 2, 4, ... , 2m, where N — 2m or N = 2m + 1, m £ N. For instance, when 
N = 4 or 5, a possible DD scheme may be visualized in terms of the following matrix, 



where each row corresponds to an even qubit and each column, supplemented with the identity operators associated 
to the odd qubits, leads to an element of the DD group, so that Q — {gj}, j = 0, . . . , \Q\ — 1 and g.j —t\® [M^+ijJa <8> 
I3 ® [M(2,j+i)]4 ® Is- Although this scheme does not scale efficiently, as the number of it pulses required to close 
a cycle grows as 4 m , it allows for the study of the effects of large control groups in DD methods with no need to 
employ excessively large systems. Numerical simulations with moderate computational resources become then viable. 
Contrary to PDD, where only one qubit is rotated at a time, the NRD sequence associated to the above Q involves 
random pulses ranging from the identity operator to m simultaneous rotations. The total number of random pulses 
leading to R simultaneous rotations is given by Qr — 3 R m\/[R\(m — R)\], where Y^r=o Qr ~ 4 m '• In large systems, 
the percentage of random pulses corresponding to a single qubit rotation, or at the other extreme, to R = m, is very 
small, decreasing with the size of the system, respectively, as 3m/4 m and (3/4) m . When m ^ 3 + 4n, the largest Qr is 
obtained for the integer R in the interval [(3m — l)/4, (3m + 3)/4], while for m = 3 + 4n, both values R = (3m— l)/4 
and R = (3m + 3)/4 lead to the two largest subsets of random pulses. 

In Fig. ^ N = 8 qubits are considered, leading to a relatively large control cycle: 256 time slots. Two situations 
favoring stochastic schemes are identified. On the left panel, the average fidelity is computed at every T n = n|£|Af. 
Even though PDD achieves first-order decoupling at these instants, the fidelity decay is substantially slower for NRD. 
This behavior persists even when the At value of PDD is shorter than that of NRD. Irrespective of the validity of the strict 
short-time condition underlying the bounds of , these findings confirm the faster convergence offered by stochastic 
methods when \Q\ is large. The fact that NRD eventually surpasses both PDD curves shows that, for sufficiently long 
times, the constraints on At for random DD may be relaxed or, equivalently, the number of control operations able to 
ensure a certain fidelity level may be smaller than in PDD. Both features may be very advantageous in realistic settings, 
given that achievable pulsing rates are finite and excessive 'kicks' might be undesirable (leading e.g. to unwanted 
heating in devices operating at dilution- refrigerator temperatures, such as quantum dots). Similar improvements are 
observed in situations where constraints on the number of pulses or control intervals make it unfeasible to close a 
complete cycle. This may be the case, for instance, when T c becomes prohibitively long. Here, no analytical fidelity 
bound for PDD exists, hence we rely exclusively on numerical simulations. When compared with the particular PDD 
sequence considered, NRD performs significantly better for most intra-cycle times t n = nAt < T Cl a result that prompts 
the search for superior deterministic sequences. Designing new stochastic methods capable of pre-filtering potentially 
good sequences simplifies this search, which would otherwise be performed over the extremely large ensemble, size of 
C(|5|™' e ')i generated by NRD. More efficient randomized protocols will be discussed in the next section. 



the DD procedure j26j 



III. CONVERGENCE IMPROVEMENT 
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FIG. 1: (color online) PDD vs. NRD based on a nested pulse sequence for Hamiltonian (3) with A = 1 and N = 8 in the logical- 
rotating frame. Left panel: average fidelity at T„ = n\Q\At, \Q\ = 4 4 . PDD: T c = 0.08J" 1 - (blue) dashed line, T c = Q.12J~ X - 
(blue) short-dashed line; NRD: (red) solid line, At — 0.12J _1 /4 4 , average over 50 realizations. Free evolution: (black) oscillating 
solid line. Right panel: average fidelity within a cycle, t n = nAt with At — 0.005J -1 for both protocols, average over 10 2 
realizations. 



IV. LONG-TIME IMPROVEMENT 



Conceptually, there are two main strategies for boosting DD performance. One rests on ways for increasing the 
averaging accuracy (hence the minimum power of At) in the effective Hamiltonian, the other on slowing down the 
accumulation of residual errors due to imperfect averaging over long times. Based on these guiding principles, we 
introduce several DD schemes and discuss their relative merits. 

In the deterministic domain, one possibility is motivated by the Carr-Purcell sequence of NMR, and consists of 
symmetrizing in time the control path of the PDD. It leads to what we call "symmetric deterministic DD" (SDD). The 
cycle becomes twice as long, but all odd order terms in Hq are also canceled. Another scheme, which generalizes NMR 
supercycle techniques 0] , corresponds to "concatenated DD" (CDD) , as recently formalized in [2]j . CDD has a temporal 
recursive structure, whose level £+ 1 of concatenation is determined by the pulse sequence Ce+i = CePiCgP2 ■ ■ ■ C^Pn, 
where Pk is the fcth pulse, Co is the inter-pulse interval and C\ denotes the generating PDD sequence. At level 1 = 1 
the concatenated sequence is also symmetric. Interestingly, however, CDD may outperform SDD even before this level 
of concatenation is actually completed, reflecting its superiority in reducing the accumulation of errors. 

In terms of randomized DD, we introduce hybrid protocols, which combine deterministic and stochastic features. 
The purpose here is to ensure good performance at short times, as typical of deterministic protocols, while at long 
times, instead of accumulating errors coherently, randomization guarantees this to happen probabilistically. The 
protocols are classified according to an inner and an outer code. The former establishes the pulse sequence in the 
interval [n|£/|Ai, (n+l)|(?|At], being associated with the control path chosen to traverse Q. Since the group is traversed 
in full, as in deterministic schemes, the inner code leads to an effective Hamiltonian H c k oc 0(Af). The outer code 
determines the additional random pulses to be applied at T n — n\Q\At. The latter are drawn from a group Q' ', which 
needs not coincide with Q. Randomization may then be associated with the choice for the inner code or the outer 
code. In the first category we have the "random path decoupling" (RPD) protocol, as proposed in 0. It consists of 
randomly choosing, at every T n , which control path to follow to traverse Q. Here, as in most stochastic protocols, 
logical and physical frame do not always coincide and we need to keep track of the applied control trajectory, so 
that an appropriate control operation may be used to correct frames. However, similarly to PDD, we may choose to 
fix the first group element as 1, which leads to frame coincidenece at every T n . This may be particularly useful, for 
instance, in conventional line-narrowing spectroscopic applications. We will call this alternative "pseudo-random path 
decoupling" (pRPD). To the second category belongs the embedded scheme (EMD), inspired to [9j. The inner code is 
a fixed PDD sequence, while the bordering pulses may either be picked at random from Q or from a different control 
set, corresponding, for example to products of uncorrelated Pauli operators as described in 9]. Contrary to RPD, this 
protocol may suffer from non-uniform performance across the set of |(?|! possible inner paths, requiring a pre-selection 
of a good deterministic pulse sequence. The performance of both EMD and RPD is significantly improved by further 
symmetrizing the inner control path in the same manner as in SDD. Here, we will be dealing only with the "symmetric 
random path decoupling" (SRPD) protocol. 
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A. Bounds on Fidelity Decay 

Analytical upper bounds on the order of the fidelity decay, 1 — E{F e (T)}, may give an insight on what to expect 
from the above protocols. Generalizing the arguments of we nn d the following order-of-magnitude estimates 

for deterministic (top line) and stochastic (bottom line) schemes: 

PDD SDD 

T 2 (\g\At) 2 n 4 T 2 (\g\At) 4 n 6 

NRD RPD/EMD SRPD 

TAtn 2 T(\g\At) 3 n 4 T(\g\At) 5 n 6 

For deterministic protocols, residual errors add coherently, which leads to a quadratic-in-time fidelity decay, 
0{{T\\H \\ 2 ) 2 ), as found in 0. Therefore, it is only the ability to cancel or reduce higher order terms in Hq that 

may induce better performance. At short times, the dominant term in each cycle of the PDD is and the bound 

is derived from the norm 

H-ffo^lb < k2t c- For SDD, = 0, and the norm of the dominant term is limited as 

ll-^o I h — k3 T 2 . In the case of CDD, the averaging accuracy depends on the level of concatenation (level 1 recovering 
the results of PDD) and on the system considered. 

Contrasted with deterministic methods, the accumulation of residual errors for random protocols is slower, as 
reflected by the linear-in-time decay of the fidelity. This may be intuitively justified as follows. Each step of NRD 
can accumulate an error amplitude up to nAt and during a time T there are T / At such intervals. Due to the 
randomization, amplitudes add up probabilistically, leading to a decay oc TAtn 2 . The reasoning is similar for the 
other protocols, but each step now corresponds to the total interval of the inner code, p|0|At, p = 1 \p = 2] for 
RPD/EMD [SRPD], leading to £{F e (T)} > 1- £)(T\G\At\ \H cS \ \\). The norm of the effective Hamiltonian is derived from 
the deterministic sequence underlying the stochastic protocol, the worst path being considered for RPD/SRPD. Hence 
H-Hefflb ~ I \Ho 1 12 and here again H-ffolh becomes the main factor differentiating the protocols. We therefore expect a 
significant better performance for SRPD (whose norm comes from SDD) than for RPD/EMD (whose norm is that of PDD). 

Merging together features of deterministic methods and pulse randomization, as in hybrid DD schemes, suggests 
the possibility of suppressing errors more effectively at both short and long interaction times. However, since the 
above bounds apply only at short times, numerical analysis becomes necessary. 

B. Numerical Results 

For the model of Eq. |0J), first order decoupling can be achieved through a very simple system-size- 
independent scheme. It consists of alternating two rotations around perpendicular axes, one acting on all 
odd qubits, the other on the even ones, such that the cycle is closed after 4 collective pulses. g = 
{1, Z\Z^ . . . Zjv-i, Zil2^3^4 • • • Zn-{Yn, Y2Y4 • ■ ■ Yn} is a possible DD group realization for even N. 

A quantitative comparison is presented on the left panel of Fig. [21 Note that the (inner) sequences characterizing 
SDD, SRPD and CDD are not necessarily completed at the instants of data acquisition, T n = AnAt. We verified that the 
outcomes of pRPD and the EMD based on a single group g are very similar, while RPD and the EMD based on random Pauli 
operators perform closely at intermediate times. RPD becomes the best of the four protocols at long times. Results are 
shown only for pRPD. As expected, random protocols surpass deterministic schemes at long times, the crossing being 
evident between protocols that have equivalent performance at short times. We note that NRD meets PDD already 
at very small values of ((F e )), since the group is now small. In addition, CDD is remarkably outperformed by the 
relatively simple SRPD method, which can be understood by re-examining the analytical bounds. For this particular 
system, CDD at level 2 achieves only second order decoupling, so that Hq ' ^ and the bound of SDD is recovered. 
It turns out that due to the reducibility of g, increasing the level of concatenation does not improve the protocol 
performance. We find that Hq contains terms such as XjXk, YjYk, and ZjZk, (j,k - odd), which are unaffected by 
the pulses drawn from g. As a consequence, performance improvement is saturated and the coherent accumulation 
of residual errors soon deteriorates the results obtained with CDD - the method is eventually outperformed by SRPD. 

Given the protocols above, a way to guarantee the best performance through the whole time axis consists of 
interpolating the CDD scheme at short times with SRPD at long times. This is illustrated in the right panel of Fig. [21 
CDD is used until the third level of concatenation is reached, at T = A 3 At, where we then switch to the SRPD sequence. 
Note that if the applied control history is recorded by means of an appropriate classical register, a randomly generated, 
optimized deterministic pulse sequence may be obtained in this way upon de-randomizing the protocol at the end. 

The identification of an efficient pulse sequence is strongly dependent on the time interval considered. Building 
on the fact that deterministic protocols perform better at short times, while stochastic schemes become superior at 
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FIG. 2: (color online) Deterministic vs. random DD based on a 4-pulse sequence for Hamiltonian JIJ with N = 8 in the 
logical-rotating frame. Average fidelity over 10 2 control realizations, at T n — n\Q\At. Left panel: A = 1, At = 0.1J -1 . Free 
evolution: (black) dashed line at the bottom. Right panel: A = 5, At = 0.05 J -1 . 



long times, it might be beneficial to exploit closed-loop strategies to determine when to switch from one to the other. 
While the process of monitoring a quantum system in real time with the purpose of controlling its dynamics has 
been proposed to a variety of systems, ranging from cavity QED |28| to nanomechanical systems p9j| . experimental 
implementations have been reported only recently |30| and the capabilities are still limited. However, assuming the 
possibility of monitoring online the system considered here, we could, for instance, decide when to switch from CDD 
to SRPD once a certain a priori stipulated value of ((F e )) is reached. 

An additional advantage of randomization, as discussed in Ref. |12(, is related to time- varying systems. In this 
case, randomized protocols lead to more robust performance, since they are usually more protected against adversarial 
situations where a pre-established control action may be inhibited by the system fluctuations. Here also, the possibility 
of real time feedback might be useful, allowing, for instance, to better adjust the pulse sequence based on the system 
parameters variations. 



V. CONCLUSIONS 



We have reinforced the advantages of randomization in terms of faster convergence and long-time stabilization, by 
comparing the performance of various decoupling schemes in refocusing the evolution of a chain of nearest-neighbor- 
interacting qubits. Our analysis also indicates the promising role of stochasticity in the search of optimized pulse 
sequences. While preliminary results indicate that the main conclusions remain unchanged when pulse imperfections 
are considered, further analysis is needed in this direction. It is also our hope that these findings will prompt 
experimental verifications in available control devices. 
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